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We discuss a simple approach to measure the weak mixing angles (fis and (/>d of the Bs and Bd 
systems in the presence of New Physics. We present a new expression that allows one to measure 
directly if New Physics contributes significantly to the mixing only. We apply the method to 
specific penguin-mediated B — >■ PP, B — > PV and B — VV modes. We provide a very stringent 
and simple bound on the direct CP asymmetries of all these modes, the violation of which is a signal 
of New Physics in decay. Within the same theoretical framework, an updated prediction for the 
branching ratio of Bs —>■ K^* is presented, which can be compared with a recent LHCb analysis. 



In our quest for New Physics (NP) in flavour processes, 
the weak angles (f>M involved in the mixing of Bd {(pd) 
and Bs {4>s) mesons have been studied with a great at- 
tention recently, since they have shown interesting dis- 
crepancies with respect to the Standard Model [l|, ■ On 
the one hand, there is a tension in the SM fit between the 
measurement of (bd and the branching ratio oi B ^ tv 
(see also Refs. 0,131 where the discrepancies with respect 
to other inputs are discussed, e.g., inclusive or exclusive 
semileptonic decays, or mixing-dependent observables) . 
On the other hand, the most recent angular analyses of 
Bs J/'Hij) by CDF, D0 and LHCb [H-Q seem to reduce 
substantially the space for New Physics in (f>s- Finally, 
there still remains an important discrepancy with the 
Standard Model through the dimnon asymmetry mea- 
sured at D0 

In this context, it is useful to devise alternative meth- 
ods to extract the mixing angle of both Bd {4>d) and 
Bs {(j)s) systems from experiment. In this paper we will 
take the point of view of assuming that New Physics con- 
tributes significantly only to the mixing phases whereas 
the formalism of the CKM matrix can still be used to 
analyse decay amplitudes. Under this assumption we 
present a clean method to extract the NP contribution to 
the weak mixing phase 0^,^ from selected Bd^s ^ PP, VP 
and VV decays. 

Some time ago [Toj . we proposed a test on the value 
of sin (ps in the SM by measuring certain longitudinal 
branching ratios and CP asymmetries of Bd.s mesons de- 
caying into vectors, for modes mediated by penguin dia- 
grams - we focused on the potential of the Bs — > K'^*K^* 
decay at that time. The method has advantages both 
from the theoretical and experimental points of view. On 
the theoretical side, it reduces the required theoretical in- 
put mainly to a single hadronic quantity A, defined as 
the difference between "tree" and "penguin" contribu- 
tions. More precisely, for a Bq meson decaying through 
a b ^ q penguin-mediated process, the decomposition 



A = A{Bq ^ Ml A/a) = Xi'^T + A^^'P , 



(1) 



that are multiplied by aI'^'' and xi^^ respectively, i.e. 



A=T - P . 



(2) 



with the CKM factors A 



can be used to 



define A as the difference of the hadronic matrix elements 



In the case of penguin-mediated decays, the evalua- 
tion of this quantity using QCD factorisation (QCDF) 
is expected to be particularly robust, as it was built 
specifically to cancel the infrared divergences coming 
from spectator-quark and annihilation contributions up 
to next-to- leading order [23| [HI. On the experimental 
side, this test of weak mixing angles within the Stan- 
dard Model can be performed simply by measuring a 
CP-averaged branching ratio and an untagged rate. This 
allows us to avoid tagging, although a time-dependent 
analysis is still required. 

In the present article, we extend this approach beyond 
the simple SM test on the weak mixing angles (jfs and 
(f>d presented in Ref. [l0| . Indeed, we show that we can 
also use penguin-mediated decays to measure the size 
of NP contributions to neutral-meson mixing, assuming 
that its contribution to AB = 1 decays is negligible. We 
also update our predictions for the branching ratio of 
Bs — > K*'^K*'~' which has been recently measured. The 
paper is organized in the following way. In Section [1] we 
present a formula to pin down the NP contribution to 4>s 
and (l)d, obtained from b ^ d, s transitions, under the as- 
sumption that NP provides significant contributions only 
to the meson- mixing phases. In Section [2] we discuss the 
main theoretical input in this method. In Section [3] we 
consider three different examples of the method, corre- 
sponding to B — >■ PP, PV and VV decays. We update 
our prediction in Ref. [3 for BRi°'^g(B^ x*°K*°) 
in Section and compare it with its recent measure- 
ment by the LHCb collaboration [l^ . In Appendix |X] 
we provide a dictionary between the experimental mea- 
surements [l^ [13] and the theoretical quantities [l^l de- 
fined for longitudinal observables m B ^ VV decays. In 
Appendix |B1 we discuss differences in the determination 
of the branching ratio for neutral mesons for tagged or 
untagged analyses, highlighting the role played by the 
width difference AF (a problem particularly relevant for 
Bs mesons). 
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1. FORMULAE FOR NP MIXING ANGLES 



Using the unitarity relation xi!'^ + Ac''' + x'f^ = 0, we 
can write Eq. ([T]) in terms of Ac'^' and Aj'' 

A = A{Bq M1M2) = -A^"'' T - A(«' A . 

The weak phase in At'' is the angle /3g, defined as 



(3) 



Pq = arg 



VtbVt 



tq 



VcbV* 



whereas Ai'^' is real to a very good approximation for 
both q = d, s. Following the definitions in Appendix El 
we introduce the observables: a branching ratio BR = 
5ps(|Ap + |Ap)/2 (where gps is the phase space factor) 
and three CP asymmetries A^ii, Amix and Aaf 



Idir 



A 



Ar 



Re(e-^'^«A*A) 



Im(e-^'^«AM) 

|A|2 + |i|2 



(5) 



defined in terms of A and its CP-conjugate A, as well as 
the Bq meson- mixing phase (pq [rjf is the CP parity of 
the final state). 

These four observables can be written, using Eq. ([3]), in 

terms of A^''] , T, A and (pQ. One can then eliminate the 
hadronic parameter T to obtain a relationship between 
these observables 

2gp,|AnA(«)|2sin2/3, (6) 

= BR{1 - 1-jf sin $QqAnix + Vf cos ^Qg^Ar) , 



with $Qg defined as 



(7) 



being (/)q the new physics contribution to the mixing 
angle of the Bq system: (f)Q = 2/3q + (f>Q^ ■ In deriving 
Eq. ^ we have assumed that New Physics could alter 
only the mixing phase of the neutral Bq meson, but not 
the CKM matrix elements involved in the decay process. 
Eq. ([6]) is a gen eralization of similar formulae developed 
in Refs. [lol in the context of the SM. This relation 
is the starting point of our analysis. 

Collecting all terms on the left-hand side in Eq. ([6|), 
and defining 



C 



2gp,|Acpsin'^,|A|' 



BR 



(8) 



we can solve Eq. ([6|) for $Qg: 



sin$Qq = Z ?//A,nix ± - ^2 77/AAr 

cos$Qq = -z rif A Ar ± Vl^z^ ?7/A,„i 



(9) 



Here, z = (1 - C)/^l - A^.^, A 
and ^Ar 



Anix/ \/l ~ ^dir 



A 



Ar 



1 — ^^ij., with the relation: 



42 



(10) 



There is a two- fold ambiguity in Eq. In practice, we 
will see that z ~ 1 (or equivalently C ^ 1), so that the 
two solutions are very close and they can be considered 
as one single solution within current theoretical and ex- 
perimental uncertainties. Notice also that z^ < 1, which 
provides a very strong and interesting constraint on the 
size of Adir: 



< C(2-C). 



(11) 



In general, a pair of values for BR and A imply an upper 
bound on \Ad[^\ [l^, HBl which is much tighter than the 
experimental value for Adir. A violation of this bound 
would be an indication of NP in decay. In the follow- 
ing we will take A^ir = ± V2C — in the numerical 
analyses. For example, in the case of i? — cjiKs we have 
C = (8.9 ± 7.7) X 10-^ which implies |Adir| < 0.019. 

If there is no New Physics and if the process consid- 
ered is such that Q = q, one recovers sin^g, — > and 
cos$Qq 1, as can be checked from Eq. ([7]). In this 
case, Eq. (jG]) reduces to the simple relation: 



(12) 



This compact equation can be easily rewritten as a SM 
test of the angles of the unitarity triangle and corre- 
sponds to Eqs. (33)-(36) of Ref. 0. 



2. THEORETICAL INPUT 

The theoretical input in this formalism is limited to 
the quantity A defined in Eq. As discussed in 

Refs. [13, [U, this quantity can be computed safely 
within QCD factorisation for penguin mediated decays 
because of the cancellation of long-distance contribu- 
tions. As a consequence of this cancellation, only penguin 
contractions contribute to A, as can be seen by inspec- 
tion of the formulae in Ref. ^| . The general form of A 
for a ^ M1M2 decay is flOlfTlj: 

(13) 

where Mi picks up the spectator quark of the Bq meson. 
The normalisation involves the M2 decay constant and a 
B — > Ml form factor: 



mljM,A''^'^''^{0) . (14) 



Ci is the relevant Wilson coefficient of the effective 
Hamiltonian at a scale of order mi,. Gm^ = Gm2 i 
Ty^'^GM^ is the penguin function with a w or a c quark 
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Channel 


|A| (10"^GgV) 


C X BR 


r>d — ? 




yZ\j.o zt zi.yj ■ lU 


JDs — r 




yL.AL zt U.oyj ■ iU 


R . /V"A 
_Dd J\(p 




l^U. / 4 zt 0.04 j ■ iU 


jDd — r xv y \ 


{ O on _|_ A QO\ 


1^14. ( zt IZ.IJ ■ lU 


R . /v^* 

£jd — r l\ l\ 


fn A^ -\- f\ ^^n^ 

l^U.'il zt U.DU ) 


/'A /17 _1_ 1 QO\ 1A~^ 
/ zt i.oo J ■ iU 


Lis — ^ 


(0 ^a -\_ f\ Qc\\ 
l^z. 10 zt u.oy ) 


l^U.OU zt U.OUj ■ iU 


Bs K*K 


(0.36 zb 0.53) 


(0.02 ±0.05) ■ 10"^ 




(1.85 zb 0.93) 


(9.37 ±9.53) - 10"^ 


Bs K*K* 


(1.62 zb 0.81) 


(0.33 ±0.33) - 10"^ 


Bd^K^cj^ 


(1.92 zb 1.03) 


(0.49 ±0.53) - 10"^ 


Bs (t>K* 


(1.87 zb 0.94) 


(8.80 ± 8.96) - 10"'^ 


Bs 00 


(3.86 zb 2.09) 


(0.92 ± 1.00) ■ 10"^ 



TABLE I: Values of A for the various decays of interest, and 
the corresponding values for C x BR, defined in Eq.®. In 
the case of two vector mesons these numbers should be un- 
derstood as referring to longitudinal polarisations (the third 
column means C X _Bi?'°"s in this case). 



running in the loop and the sign depends on whether Mi 
is a pseudoscalar or a vector meson (see [H, |T3] for the 
precise definition of the quantities entering A). For iden- 
tical particles in the final state (for instance Bs — )■ (t><f) 
an extra factor of 2 must be inckided in A. 

In the case in which both Mi and M2 are vector 
mesons, there is a different A for each of the three 
transversity amplitudes. As discussed in Ref. [l^l, only 
longitudinal amplitudes can be handled safely within 
QCD factorisation, so we shall focus on longitudinal 
polarisations for vector-vector modes, in the spirit of 
Ref. (see also Ref. ^ for a discussion of trans- 
verse amplitudes in QCDF). This requires the definition 
of suitable longitudinal observables, together with a pro- 
cedure to extract them from the observables that are ac- 
tually measured, as discussed in more detail in the ap- 
pendices of the present article. 

In Table U we present the values of A for Bq — > XY 
penguin-mediated decay channels, with X,Y — P, Vl 
(where Vl stands for longitudinally polarised vectors), 
using the same hadronic inputs as in Ref. [lOj . 



3. NEUTRAL-MESON MIXING ANGLES FROM 

Bq PP,PV AND VlVl 

In this section we illustrate the determination of the 
NP contributions to both weak mixing angles (j)d and 0s 
using Eqs. ([5]) and (jH]). We consider channels of three 
different types: Bs K"K° {B PP), Bd (j^Ks 
(B PV), and Bd ^ c/^K* {B VlVl). We will see 
that in the case of Bd — > (/>i^s and Bd — J> 4>K* , all the 
required observables are already measured, which allows 
us to perform a complete analysis to extract . There 
is no experimental data for Bs — !• A'°A'", so we will only 



be able to present a case study, which can be exploited 
immediately as soon as data is available. 

Besides the theory input A and the measured branch- 
ing ratio, we need additional inputs for the SM CKM 
parameters. Since we are assuming that meson mix- 
ing could be affected by NP, wc must use values of the 
CKM elements extracted using modes insensitive to mix- 
ing. The fit in Table 11 of Ref. [l| satisfies this require- 
ment, as the mixing-related observables included in this 
fit are used essentially to determine the size of NP in me- 
son mixing and have a very small impact on the CKM 
parameter themselves, which are thus determined from 
tree-dominated quantities. Indeed we checked that a fit 
limited to the tree dominated inputs of Ref. [l| yields 
almost identical results for the CKM parameters. From 
Ref. [l| we obtain, upon symmetrizing errors: 





= (4.1 ±0.3) • 10~^ 




= (9.5 ±0.6) • lO^-* 




= (9.2 ±0.2) • 10"^ 




= 0.041 ±0.001, 


Pd 


= (26.2 ±2.1)°, 


f3s 


= -(1.26 ±0.07)°, 
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= (69.9 ±4.4)°, 



which will be used in the following studies. 

3.1. First example: Bd — ^ 0Jfs 

We first consider Bd — > (/>A's. Using the formulae 
in Eq. 1^, together with A and C in Table U and 
the experimental numbers BR = (8.3 ± 1.0) • 10~^ and 
Anix = — 0.56 ± 0.16 [1^ (the product of intrinsic pari- 
ties is r/0Xs = ~1): we find the two solutions for the NP 
contribution to the mixing angle: 

(t^T = (-0.38 ± 0.21) V (-0.35 ± 0.21) rad , (16) 

from which wc can give an averaged result: 

< aver = -0.36 ±0.22 rad. (17) 

This result is (marginally) compatible with the SM value 
^ = 0. However the error on (fy^^ is almost completely 
dominated by the experimental uncertainty in Amix^ any 
improvement on the latter would impact our knowledge 
of (f)^^ and its agreement with SM expectations. 

In Figure [T] we show the regions in the Amix-^Ar 
plane corresponding to different values of <p^^ j the NP 
contribution to the mixing angle. The easiest way to 
understand this plot is through the two constraints 
on ^niix and A/y^r'- Eq. (jS]) is a linear equation at 
fixed (j)Q^, whereas A^^^ + Aj^-p < 1 yields a radial 
constraint. For a fixed 0q^, the solution for ^mix and 
is thus given by the intersection of a line with 
the unit circle. If this line is tangent, the solution is 
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FIG. 1: Ajnix vs. ^Ar, for several values of the NP mixing 
angle 4>Y for Ba 4>Ks, where BR = (8.3 ± 1.1) ■ 10"^ 
A = (2.32 ± 1.00) ■ lO-'^GeV and gps = 8.4 ■ 10^ GeV"^ 
The box indicates the experimental value for the asymmetry 
Amix ~ —0.56 ± 0.16 and the range obtained for Aaf using 
^dir = 0±0.019 inf = -1). 



limited to one point, and our determination is perfectly 
precise. The opposite situation occurs if the line goes 
through the origin. It turns out that the distance of 
this line from the origin is given exactly by 1 — C, so C 
measures the precision with which we can determine (j)Q^ : 

The smaller the value of C , the more precisely we can 
pin down the value of (pQ^ . 

The SM solution is shown explicitly in Fig. [1] and the 
gray box indicates the current experimental situation. 
One can see there is marginal agreement, in line with 
the result in Eq. ([T7| . We stress that these results in- 
clude all known hadronic uncertainties, and the errors 
are relatively small because the chosen theoretical input 
is robust. A more precise value for the angle Pd would 
achieve a substantial reduction of the regions of fixed 
in Fig. m 

It should be mentioned that, since we have no way 
of determining the sign of Aaf, there is a second box 
in Fig. [T] that has not been drawn, corresponding to 
T]fA^T ~ -1-0.8. We have discarded this possibility be- 
cause it leads to a very large New Physics mixing angle. 



3.2. Second example: Bd ^/T* 

An angular analysis of the decay Bd — > <I>K* is avail- 
able from both Babar and Belle collaborations , 
with an additional time-dependent analysis for the for- 



FIG. 2: Amix vs. Aaf, for several values of the NP mixing 
angle (^"^for Ba (t)K* , where = (4.7 ± 0.4) ■ 10~^ 

A = (1.92 ± 1.03) • lO"'^ GeV and gps = 8.2 ■ 10^ GeV"^ The 
box indicates the experimental values for the asymmetries 
given in Eq. ^ {rff = +1). 

mer experiment. The averaged results read [Tsl. [l9| 

BR = (9.8±0.7) • 10-^ Acp = 0.01 ± 0.05 , 
/l = 0.48 ±0.03, A^cp = 0.04 ±0.06, 

A(/.o = 0.28 ±0.42, A<5o = 0.27 ±0.16. 

(18) 

Using equations (jA9p - (jA12p in Appendix 1X1 (and noting 
that ri(f,K* = +l)j we obtain the longitudinal observables: 

Al°^^ = -0.05 ±0.08, 

A^°'-^ = 0.96 ±0.25, (19) 

A^^^ = ±(0.27 ± 0.86) . 

Compared to the case of Bd — > 4>Ks^ the uncertainty on 
Al^^ is so large that one cannot distinguish, among the 
two solutions, which one is disfavored by too large a value 
of the NP mixing angle. For illustration we focus on the 
negative solution of Al^. 

For this decay, Eq. (fTTj) yields the following range for 
the direct asymmetry: 

^k>ng ^ Q _^ Q ^ ^20) 

which is more precise than the experimental value in 
Eq. (jl9p . and which is used in the following. Equation 
([9]) then yields the two solutions: 

= (0.31 ± 0.90) V (0.34 ± 0.90) rad , (21) 

which can be averaged as 

'/'raver = 0.33 ±0.90 rad. (22) 
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This result is compatible with the SM (0^^ = 0), and 
also with the result obtained from Bd 4>Ks (Eq. (|17p ). 
within large uncertainties. 

In Fig. [21 we show the regions in the Amix-^Ar plane 
corresponding to different values of the NP contribution 
to the mixing angle. The SM solution is shown explic- 
itly, and the gray box indicates the current experimental 
situation. From this plot we see that a more accurate 
knowledge of ^I™^ or A'^p^ would be very useful to test 
the SM hypothesis. 



3.3. Third example: Bs K°K° 

There is currently no information on this mode. How- 
ever, the U-spin related decay Bd K^K^ has the mea- 
sured branching ratio [l9} : 

BR{Bd -> K°K°) = (0.96 ± 0.20) • 10"^ (23) 

As explained in Ref. if there is no NP affecting the 
decay Bd — > K^K^, one can predict the SM branching 
ratio for the decay Bg — > K^K^ quite precisely. This 
prediction is not affected by a possible NP contribution to 
4>d (see Section 3 of Ref. and since we are assuming 

that there is no NP affecting Ai? = 1 processes, it is a 
good option to take this prediction as an input for our 
analysis 

BR{Bs K°K°) = (18.2 ± 7.3) • 10"^ . (24) 
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FIG. 3: Amix vs. Aaf, for several values of the NP mixing 
angle cj>'^^ for Bs — ^ K'^K'^, where we take the estimate BR = 
(18.2 ± 7.3) • 10"'' (see the text). In addition, we used A = 
(3.05 ± 1.11) • 10"'^ GeV and gps = 8.03 • 10^ GeV-^ 

The three plots in Figs. 11121 and [51 can be reversed and 
reinterpreted as predictions for the ^mix and Aaf asym- 
metries if the range for cj)^ ^ is extracted from other pro- 
cesses, always under the assumption that there is no sig- 
nificant New Physics contribution to Ai? = 1 processes. 



The present analysis of this decay should be considered 
as a case study, waiting for data to be available. 

Without further experimental information on CP 
asymmetries, we cannot use Eq. © to extract (j)^^ , the 
NP contribution to the mixing angle. But it is still pos- 
sible to determine the regions in the Anux'-^Ar plane 
that correspond to each value of 0^^, similarly to Figs.[Tl 
and [21 Notice that in the present case we have simply 
= 0NP^ gjj_^(,g Q ^ g i^cre. The case Bg ^ K^K^ 
is shown in Figure [31 Whenever the branching ratio for 
this mode is measured, this plot can be remade, although 
we do not expect it to change in an appreciable way. A 
precise measurement of the time-dependent CP asymme- 
try will then provide an accurate determination of . 

One may wonder why we did not consider the decay 
Bd — >■ K^K^ at a first stage, since its branching ratio is 
known. However, the prospects for this particular decay 
are not very alluring. The mixing angle involved in the 
expression for C (Eq. ([5])) is jSd, and the branching ratio 
is not particularly large, so that the value of C is small 
but not tiny (C ~ 0.03). As mentioned above, the value 
of C determines the potential accuracy of our procedure, 
which turns out to be rather poor in the present case. 
This is true in general for all & — >■ d transitions in our list 
(see Table HI; however, the situation may change if the 
branching ratio of some of the modes is found to be large 
enough, as it might be for Bd — )• K*"K*°. 



4. THE BRANCHING RATIOS BR{Bd,s ^ K*°K*°) 
WITHIN THE SM 



In this section, we present a prediction for 
^ K*°K*°) in the SM using our approach. 
This prediction can be easily turned, as shown later on, 
into a prediction for the total branching ratio, once /l, 
Ac p and A^p are measured (see Appendix [X] for de- 
tails). This section is an update of results presented in 
Ref. [3. 

The two decays Bq — )• Ar*°A^*°, with q — d,s, are 
related by U-spin symmetry. The symmetry-breaking 
effects can be separated into a factorizable and a non- 
factorizable part. This translates into the following rela- 
tions involving the hadronic parameters Tq and Pq 



Ps = f Pdil + Sp) , Tg^f Tdil + St) 



with the factorisablc factor / 

f 



K'K* 



(25) 



(26) 



and 5t,p account for the non-factorizable symmetry 
breaking effects. There is no theoretically clean way 
to compute these quantities. Here we use QCD- 
factorisation to estimate an upper bound on this cor- 
rections. We get \5p\ < 0.09; and keep its phase as a free 
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5 • 10-'^ 


5 ■ 10-"^ 


5 • 10-5 


ODMV 


16.05 ±4.87 


16.38 ±4.92 


16.46 ± 4.93 



TABLE II: Results for the ratio R™^ for three different val- 
ues of BR'°''<'{Bd K'°K*°). It can be seen that the de- 
pendence on this branching ratio is very mild. 



parameter (even though QCD factorisation would predict 
it to be small) . It is very easy to show that the quantity 
6t is related to Sp through 



Td St ^ PdSp . 



Observables related to the decay Bs — >■ K*'^K*° can be 
obtained if the ones for Bd x*°K*° are known. An 
important result discussed in Ref. [l^ is that the ratio 



Rsd = 



_ ^ K*°K* 



(28) 



is almost a constant if BR^°'^^{Bd K*°K*°) > 10"^, 
even in the presence of NP in mixing. Indeed, when 
the two branching ratios in Eq. ([28| are reexpressed in 
terms of Ps^d and Ag^d (as computed in Table H]), it be- 
comes clear that the ratio Rsd is completely dominated 
by the penguin contributions Ps and Pd for large enough 
BR''°^9(^Bd K*°K*°), and thus determined by / and 
Sp essentially. This behaviour is confirmed in Table |TT] 
for several values of the branching ratio. 

Our result, combining all the error sources, is 

= 16.4 ±5.2 , (29) 

which updates Table IV in Ref. • The improvement on 
the uncertainty comes in particular from Eq. (|27p which 
was not used in this reference. We have added the su- 
perscript "DMV" to distinguish our determination from 
other ones that we will describe now. 

Indeed this number can be compared with the value 
obtained within the QCDF framework taking the usual 
model for l/m;,-suppressed corrections described in 
Ref. [Ol. The corresponding predictions are for the Bg 
decay mode: BR{Bs ^ K°*K°*) = (Q.ll^^^g^) x IQ-^ 
and fhiBs) = 0.63to'29i and for the corresponding Bd 
decay mode: BR{Bd -> K°*K°*) = (0.6ioj) x IQ-^ 
and fL{Bd) = O.GOlo'jy- We write Rsd in terms of total 
branching ratios and polarisation fractions, as described 
in Appendix \K\ 



_ BR{Bs)fUBs) . 
"^'^ BR{Bd)fL{Bd)^'' 

where we introduce the correcting factor 



/c 



l + A%p{Bs)Acp{Bs) 
l + Alp{Bd)Acp{Bd) 



(30) 



(31) 



In Ref. 

Acp{Bs) 



[iTf, the predictions within QCDF arc 
± 1% a.nd Acp{B,i) ~ -13%. Since the QCDF 



predictions are dominated by longitudinal polarisation, 
we assume A^p{Bs) and AP(jp{Bd) to be of the same 
size, leading to a factor /c ranging from 0.98 to 1, and to 
the QCDF prediction 



= 13.8 ± 19.2 



(32) 



Alternatively, if data on B ^ K*<f> is used to con- 
trol annihilation rather than the usual model for these 
1 /TOfc-suppressed corrections, the predictions become 
[i3: BRiBs ^ /?o*ifO*) = (7.9l^j) x 10"^, /l(S,) = 
0.72toll., BR{Bd ^ K°*K°*) = {O.Gtoi) x IQ-^ and 

-1+0.16 



(27) fhiBd) = 0.69Ig'2o. In this case the prediction is practi- 



cally the same but with smaller errors 



pQCDF- 
^sd 



= 13.7 ± 10.5. 



(33) 



which is consistent with our result because of the large 
error bars, even though the central value is a bit low 
compared to ours. 

The most recent experimental values [l§| for BR{Bd ^■ 
/?o*^o*) ^ (1.28i°jg ± 0.11) X 10-6 and /^(Bd) = 

0.80t{j:}^ ± 0.06 are consistent with the QCDF pre- 
diction, although both are on the high range. The 
LHCb collaboration has measured the Bg mode very re- 
cently, reporting the following numbers [l2|: BR{Bs — > 
K°*K°*) = (2.81 ± 0.46 ± 0.45 ± 0.34) x 10"^ and 
/^(B,) = 0.31 ±0.12 ±0.04. 

Following the discussion in Appendix [B] we can ex- 
press the ratio Rsd in terms of experimentally measured 
quantities: 



Rsd 



BRLHCb{Bs)fL,LHCb{Bs) 1 - 

BRb — fact {Bd).fL,B-fact{Bd) 1 +ycos( 



(34) 



where we define y = ATs/{2Ts)- In both cases, it was 
assumed in the experimental analysis that there is no 
CP- violation in decay, and therefore we do not include 
the corrections associated with the direct asymmetries. 
We take y = 0.046 ± 0.027 [3, and the SM value of 0^. 
These results imply a quite low value for Rsd- 



R 



exp—I 



.1±3.3) X 



hjBs) 
0.31 



.1±4.7. (35) 



The measured value of fL{Bs) is unexpectedly low with 
respect to the polarisation fraction for the C/-spin re- 
lated channel: /l(B<j) ~ 0.8. If /l(B,s) ^ fUBd), we 
would obtain a ratio Rsd ~ 21 in better agreement with 
Eq. (1211), although /l(Ps) ^ 0.6 would fit better. 

An alternative determination of this branching ratio 
using data from B^ — >■ Dfn^ yields the result 
BR{Bs K°*K°*) = (2.64 ± 0.61 ± 0.42) x lO''^ in 
fair agreement with the previous result: 



Rlf-'' = (7.6 ± 3.2) X 



'fL{Bs)\ _ 



0.31 



= 7.6 ±4.5. (36) 



A similar exercise with fL{Bs) — 0.8 yields Rsd — 19, 
once again in better agreement with our expectations. 
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Finally, one can invert these relations and get a pre- 
diction for the total BR{Bs — > K'^*K^*), as measured at 
a hadronic machine (thereby the subscript 'LHCb'): 

BR{B, ^ K'*K'*rLlibb = 

fL,B~fact{Bd) 1+2/ cos 0s 



(5.7 ±2.5) • 10~^ X 



fL.LHCb{Bs) 1 - 
0.31 



fdBs) 



(5.7 ±3.4) • 10' 



(37) 



which can be trivially updated if the longitudinal polar- 
isation of the Bs mode changes. 

We remind the reader that, as opposed to the rest of 
the paper, all the predictions in this section have been 
obtained within the SM. A detailed analysis of the impact 
of NP on these observables is certainly worthwhile, but 
lies beyond the scope of the present article. 



CONCLUSIONS 

We have shown how to exploit B PP, PV and 
VV penguin-mediated decays in order to obtain the New 
Physics contribution to the weak mixing phase of both 
Bd and Bs systems, under the assumption of no signifi- 
cant New Physics contribution to the decay amplitudes. 
The main theoretical input consists of the infrared-safe 
quantity A, which can be evaluated within QCD factori- 
sation with a good theoretical control. It is interesting 
that the experimental knowledge of the branching ratio 
by itself yields a stringent bound on the direct CP asym- 
metry, Eq. ((TT|) . As a consequence, the knowledge of BR 
and Amix is sufficient to determine the New Physics con- 
tribution to the weak mixing angles. Different examples 
are discussed and the corresponding values for (pQ^ pre- 
sented; however, more precise data is required before a 
clear conclusion can be reached. 

These results should be compared with the one ob- 
tained from charmonium modes, where one is also sensi- 
tive to the NP mixing angle. Time-dependent analyses of 
B J/fA's give the average (jjf = (21.2±0.9)°. 

Using the value quoted in Eq. ([TSl) for the SM contribu- 
tion /3d leads to the estimate 



/NP 



(cc) = -0.09 ± 0.04 rad 



(38) 



In Ref. [l| , the analysis of NP contributions to both Bd 
and Bs mixings (so-called scenario I) including the ex- 
perimental information available at that time led to the 
following value for the NP mixing phase in the Bd sector 



-0.22t"o:°^ rad 



(39) 



The main ingredient for this value is the discrepancy be- 
tween 0^"^ and BR{B — >■ tv) in the SM, and thus it will 
only be marginally affected by the CDF, D0 and LHCb 



(cc) 
Seen. I [1 

Ba 4>K 



9d 



FIG. 4: Comparison between the NP mixing angles obtained 
from cc, from Bd — >■ ^^s, from Bd — >■ ^Ji'*, and from the fit 
in Scenario I of Ref. The line cp^^ — corresponds to the 
Standard Model. 



updates concerning the Bg sector. The comparison of 
these numbers with the ones obtained in Section [3] is 
shown in Fig. 2] 

Finally, an updated result for the longitudinal observ- 
ables of the decay mode Bg — ?> A'°*Ar°* has been pre- 
sented. For the ratio of Bg to Bd longitudinal branching 
ratios, defined in Eq. ([28| . we find: 



16.4 ±5.2. 



(40) 



This result can be compared with recent experimental 
analyses, and with similar theoretical predictions ob- 
tained with alternative approaches: 

1. Our result is compatible with the QCDF predic- 
tion (sec Eq. ([32])), although our error is about four 
times smaller. If the QCDF is supplemented with 
Bd — >■ K*(j) data to fix 1/mh-suppressed corrections 
of phenomenological relevance, the error decreases, 
but still is twice as large as ours. 

2. Our result is compatible with the recent LHCb 
measurement, Eq (j35p . although the experimental 
number seems rather low. We tentatively identify 
this anomaly with the surprisingly low value for the 
longitudinal polarisation fraction /^(i?^) — 0.3: a 
value of ~ 0.8, similar to the one measured 
for the J7-spin related Bd mode, would lead to a 
Rsd — 20, more in agreement with our predic- 
tion. If the LHCb measurement of /^{Bs) is con- 
firmed, its interpretation would constitute a theo- 
retical challenge, as it would require a type of NP 
that violates the flavour symmetry relating Bd and 
Bs decays very significantly. 

We hope that these results will trigger more precise 
experimental analyses of penguin-mediated decays both 
in Bd and Bg sectors, considering the potential of these 
channels to identify New Physics in meson mixing. 
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Appendix A: Dictionary for longitudinal observables 
in B — > VV decays 



In the case of two-body B decays into vector mesons 
Vi and V2 , there are three different configurations of he- 
licity available h = +1,0,-1. This is described by three 
transversity amplitudes ^o,-L,|| corresponding to linearly 
polarised states, with definite CP properties (Ao,|| cor- 
respond to a CP-parity 770,11 = 77^1^2; whereas Ai has 
77 j_ = —riViVV2)- A-S discussed extensively in Ref. [l7t[l8|. 
the transversal amplitudes are suppressed by pow- 

ers of l/ms compared to the longitudinal amplitude Aq, 
leading to an expected longitudinal polarisation close 
to 1. In turn, this implies that predictions based on the 
heavy-quark limit, such as those from QCD factorisa- 
tion or SCET, are much more reliable for longitunal ob- 
servables than for parallel and perpendicular ones, where 
only rough estimates can be derived. 

This theoretical situation explains why we considered 
only longitudinal observables for VV modes in Ref. [l^ 
as well as in the present article. We defined longitudinal 
observables as the CP observables constructed by con- 
sidering as final CP eigenstate the two-meson state with 
longitudinal polarisations, i.e., the CP-averaged branch- 
ing ratio and the CP-asymmetries 



BR 



lone 



long 



A 



/(long , .llong 



dps' 

\Ao\ 



\Ao\' + \Ao\' 



2 



I An 



-2770 



l^ol 

-"t"^A*A, 



0^0 



(Al) 
(A2) 

(A3) 



with the phase space factor relating an amplitude to the 
corresponding branching ratio 



spect to the usual angular variables 61,62,1^ [T7| ) 
B„^ViV2 9 



d cos 6*1 rf cos 92dil} 



IT 

^ ^- sin^ 01 SIT? { 



-[|ao| cos 9 1 cos 62 



,1 



I cos ij} 



sm Ui sm 
, 1 



'2 sin -0 
sin 26*1 sin 202 cos i> 



2V2 

Im[aoa^] — -= sin 26i sin 2^2 sin -0 
2v 2 

Im[a|'|aj^]— sin^ 9i sir? 62 sin 20] . 



(A5) 



If the initial meson is charged, or if we neglect the differ- 
ence of lifetimes between neutral mesons (AF = 0), the 
coefficients and the transversity amplitudes Ai can be 
identified up to a global normalisation factor, and the 
quantities (|Aip - (|A3p are easy to obtain from an angular 
analysis of the decay, as explained in the following. 

For the moment, the only available detailed time- 
dependent analyses for VV modes of interest to us con- 
cerns Bd — >■ (j^K* [l3, US] , providing CP-averaged param- 
eters and CP-asymmetries for the 3 transversity ampli- 
tudes (these analyses for the B^ meson set AFb^ = 0). 
Using a — superscript for B^ observables and a -I- super- 
script for Bd observables, one has values for the branching 
ratios 



BR- = 



T 



Ftotal 

F 

Ttotal 



= 9ps \A> 



X 
X 



\Ax\'.. 



and for the polarisation fractions 



ft- 



We have therefore the relationships 



(A6) 



(A7) 



1 1 



BR = -- 
2F 



total 



(f + F) 



fL = ^{.fl + .fL) 



Acp = 

A° - 
j\cp — 



F-F 

fTf 

ft 



fl 



(A8) 



ft + fl 



and we can easily define the longitudinal observables in 
terms of these observables 



dps 



(777 - mv^y 



(A4) 

Aq represents the longitudinal decay amplitude for Bq, 
Aq its CP conjugate, and (f>Q the meson-mixing angle. 

In principle, these quantities can be obtained from the 
angular analysis of the differential decay width (with re- 



Bi?'°"s = Bi? • fL ■ [1 + A°cp ■ Acp] 



/(long _ 



A" 



Acp 



14 



Acp ■ Acp 



(A9) 
(AlO) 



(A^°f)2sin(2/3 + arg(AoMo)), (All) 



AAr' = -V\^l (4D^ cos(2/3 + arg(AoMo)) , (A12) 
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where the relative phase between Aq and Aq is defined 
following the notation of Ref. [l3| 



arg(AoMo) = 2A<5o + 2A0o. 



(A13) 



In this analysis, the mixing angle (j)d was assumed to 
have its Standard Model value 2/3 = 0.75 ± 0.03 rad, 
which we used to translate the measured quantities into 
the longitudinal observables A^°'^^ and A^""^^. In the ex- 
perimental analysis of Ref. [13| , there is only sensitivity 
to sin(2/3 + arg(ylo/ylo)), which means that there is a sign 
ambiguity in ^^p^, as seen in Eq. ([T9|) . 



Appendix B: Neutral-meson observables with a 
finite width difference 



Since we consider the meson in the present arti- 
cle, we have to discuss the modifications induced by a 
finite lifetime difference, before distinguishing its impact 
on flavour-tagged and flavour-untagged analyses for the 
longitudinal observables of interest. We will see that 
we must include different 0{ATs/Ts) corrections in each 
case to connect the measured quantities and the observ- 
ables (jAip - (jA3p . As will become clear in the following, 
most of this discussion applies not only to two-body vec- 
tor modes but more generally to any decay of a neutral 
meson into a CP eigenstate. 



1. Branching ratios in presence of meson mixing 

The B-B systems can be described in terms of CP- 
conjugate flavour states \Bq) and \Bq) [Q = d,s). The 
time evolution of an isolated neutral Bq meson of a given 
flavour at t = and decaying at a later time t into a CP 
eigenstate / is given as [2l|, [22 



T{BQ{t) -^f) = Nf'- 



cosh 



Art 



-|-^dir cos(AAft) — Aaf sinh 



Art 



Ajjjix sin Ail/t 



(Bl) 



\Af\ 



T{BQ{t) ^ f) = Nf 



-^dir cos(AAft) — ^Ar sinh 



\Af\\-rt 



cosh ■ 



Art 



Art 



(B2) 



where the Hamitonian eigenvalues ij-l.h = 1^Il,h — 
iTL^n/^ define Afi = fin - t^L = AM - iAT/2. Also, 
r = (r^ -I- rL)/2 and Aj = (/|Bq). is a time- 
independent, but state-dependent normalisation factor, 
corresponding to the integration over phase space. The 
mixing ratio q/p = exp(— z^q) is assumed to be a pure 
phase in the present article (as suggested by the very 
small values of the flavour specific asymmetries both for 
Bd and Bs (3). 



The normalisation factor Nf comes from 

TiBQit) ^ f) = Nf\{f\BQm\ (B3) 

Let us notice that this "definition" is rather ambiguous, 
as one generally considers initial and final states that are 
asymptotic (mass eigen)states, which Bq is not. How- 
ever, one can still determine this factor by going back 
to the derivation of Fermi's golden rule. This is gen- 
erally done for an initial mass eigenstate of the unper- 
turbed Hamiltonian but can be easily adapted to the 
superposition of two mass eigenstatcs \Bi,),\Bh), pro- 
vided that the difference of energy between the two tran- 
sitions ujl — ^^H = {f\Hi\BL} — {,f\Hi\BH) is small. In 
that case, the normalisation Nf corresponds to the phase 
space available, computed for the mass eigenstates of the 
unperturbed Hamiltonian. We obtain therefore the same 
normalisation as in the case of charged i3-decays 



Nf = gpsT , 



(B4) 



where gps is given by the phase space with an incoming 
meson of mass M = {Mh + Ml)/2. 

Since there is no unambiguous definition of the CP- 
averaged branching ratio for neutral mesons as the states 
involved {Bq and Bq) are not mass eigenstates, we 
should define what we call the CP-averaged branching 
ratio for Bq decays. We opt for the simple definition, in- 
spired by the charged meson case, and corresponding to 
the value that we would obtain through a measurement 
at t = (i.e., before neutral B-mcson mixing could take 
place) 



BR 



f 



gps{\Af\' + \Afn/2. 



(B5) 



We will show later that this definition coincides exactly 
with time-integrated CP-averaged decay widths in the 
limit AT 0, but that the relationship is corrected by 
terms of order 0(Ar/r). This correction depends on the 
experimental setting, because each one has a different 
sensitivity to the time evolution of the neutral _B-meson. 
This correction factor can thus be seen as a correction for 
the temporal acceptance of the considered experiment. 

Let us now come to Bq — )■ V1V2 decays, and let us 
assume that all the direct asymmetries vanish, whereas 
^o,mix = ^11, mix = -^±,mix = ?7 shi and ^o.Ar = 
— rj cos 4>Q, where 77 = Tyy^Tyy^. As 
we can obtain the time-dependent 
5]), upon the identification 



^11, Ar ~ — ^±,Ar = 
discussed in Ref. [2l| . 
decay width from Eq 



\af\'' ^ N\Af\^e 



cosh ■ 



Art 



rjf cos (f)Q sinh ■ 



Art 



rjf sin( 



(B6) 
I sin Amt , 



Re[a;a||] ^Ar|^oP|||cos((5|| -(5o)e" 

, Art , . , Art 

cosh h rj cos (pQ smh 



-rt 



(B7) 



rj sin (pQ sin Amt 



Im[a}a_L] ^ N\Af\\A±_\e~^^ x sin(5/ cos(Amt) (B8) 

Arti 

—7/ cos 5 f cos (pQ sin( Amt) + 77 cos 5 / sin (pQ sinh — — 
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with / = 0, _L, II, and Sq and the relative strong phases 
of Aq and Ay with respect to A±. There is a common 
normaUsation factor N ~ gpsT for all amplitudes. The 
CP-conjugate expression can be obtained by the replace- 
ment (f>Q — —(f>Q and by multiplying by (—1) the two 
imaginary parts involving a± in Eq. 



2. Tagged analysis 

Flavour-tagged analyses are particularly easy to per- 
form in B-factories as they produce intricated BB 
pairs [ill, and they allow one to study separately the 
time-dependence of B and B samples as well as to extract 
the modulus and relative phases of the three transvcrsity 
amplitudes through an angular analysis. The probabil- 
ity of the process depends on the decay times itag and 
tcp of both mesons (the one decaying into a tagging 
state /tag and the one actually studied for CP-violation 
and decaying into fcp)- After integrating over the sum 
^tag -\- tcp, one obtains for rtag(Sg(0 ~^ /) a structure 
similar to that for untagged decays provided that the 
exponential exp(— Ft) is replaced by exp(^F|t|) and that 
t = tcp — ttag is allowed to run from — oo to -foo. Indeed, 
from Ch. 1 in Ref. HI] , we have the joint decay amplitude 
A(ttiig,tcp) which can be integrated over time to yield 



dt 



tag 



t: N, 



dtcP / [dp\A{Uag,tcp) = 2^1^ 



nag I 



dtN, 



I A 



cp- 



CP\ 



\A 



cosh - 



AFi 



-Adir cos(Amt) — Aaf sinh 



AFt 



Amu sin Ami 



(B9) 



= -Si?(BQ ^ /tag) X / dt 

tagged 

where the first factor comes from the angular integral de- 
scribing the e^e'' — > BqBq transition, and we have de- 



fined Ftaggod(SQ(t) — > fcp) as the integrand of Eq. (|B9 

In the case of a non-vanishing width difference, and 
assuming the normalisation Eq. (jB5[) . we can determine 
the CP-averaged branching ratio by considering 



BR 



•/.tagged 



/+oa 
[rtagged(SQW^./) 
-OO 



+Ftagged(BQ(i) ^ /)] = BR 



ThTi 



(BIO) 



We have therefore a correction of the time-integrated 
branching ratio with respect to the branching ratio in 
absence of mixing. This correction is due to the differ- 
ence of the widths between the two neutral states. This 
is typically a small correction: if we define y = AF/(2F), 
we can write UM 



L 



l + (0±2 

1 + (2 ± 2) • 10-3 for B. 



10- forS,, ^g^^^ 



3. Untagged analysis 

In the case of hadronic machines, such as CDF, DO and 
LHCb, we encounter a rather different situation with no 
fiavour tags available and an integral over time being 
performed. Since there is no information on the second 
_B-meson being produced, one must consider Eqs. (jBl[) 
and (|B2[) for t > Q. The CP-averaged branching ratio 
can then be determined through 



BR 



/,untag 



r 1 

j ^[TuntMit] ^ f) (B12) 

AF 



-runtag(SQ(t) ^ /)] =BR 



Li H 



I- A 



Ar 



2F 



Compared to the tagged case, we see that there is a 
further term coming from Ay Ar- The resulting correc- 
tion for the branching ratio is larger than the one for B- 
factories, since it is linear in the small quantity AF/(2F), 
and not quadratic. 

For vector-vector channels, assuming that the produc- 
tion rate for Bq and Bq is the same and that there is 
no CP-violation in decay, -Bi?ViV2,untag yields an angular 
structure identical to Eq. (|A5|) . with the replacement 



FrF 



Li H 



1 



AF 

'//^cos< 



Re[aoa||] gp^l Ao| |A|| | cos((5|| - (5o) 



FlF 



Li H 



1 



AF 

'/^COS( 



Im[a^aj.] gps\Af\\A±_\ 
F2 



FlF 



H 



Vf^cosSf smc 



,(B13) 



(B14) 



(B15) 



where the cos(j) terms are the remnants of Aaf- After 
integration over the angles, we obtain the CP-averaged 
width 



^Bq^ViV2 = 9ps 



AF 



F3 



|Aop + |A||p + |Ai|2 



+'?^cos0q(|Ao|^ + |A|||^-|A^|^) 



(B16) 



which is expressed in terms of the longitudinal and trans- 
verse CP-averaged branching ratios. 

In the case of the LHCb analysis for Bg j{*Oj^*o 
(77 = 1), an angular analysis was performed to extract the 
total branching ratio as well as the longitudinal polarisa- 
tion. The expression in Ref. for the latter quantity was 
given in the SM case with a vanishing mixing angle (j)s- 
However, the angular analysis performed does not rely 
on this assumption, so that we can just identify the coef- 
ficients with the same angular dependence in the differ- 
ential decay width, leading to: fh.LHCb = |aol/rB-j-ViV2 ■ 
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We can thus derive the relation, vahd in the absence of 
CP-violation in decay 



fL,LHCb X BRLHCb — BR 



long_ 



Li H 



1 



AT, 

2r. 



■ cos ( 



(B17) 



We see that even in the SM case where 4>s is tiny, there 
is a small contribution from AFs/Fs to the relationship 
between the polarisation measured at LHCb and BB}°^^. 
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